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ABSTRACT 



u 

^ : 

*T^' We have carried out two extremely deep surveys with SPIRE, one of the two cameras on Herschel, at 250 /im, close to the peak of the far-infrared 
q background. We have used the results to investigate the evolution of the rest-frame 250-yum luminosity function out to z = 2. We find evidence for 
5^ ■ strong evolution out to z — 1 but evidence for at most weak evolution beyond this redshift. Our results suggest that a significant part of the stars 
"^? ' and metals in the universe today were formed at z < 1.4 in spiral galaxies. 

i ^| Key words. Galaxies: evoluton - Galaxies: formation - Galaxies: high-redshift - Submillimetre: galaxies 

^H ■ 
^ 1 . Introduction significant fraction of the FIRB at wavelengths where its energy 
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density is at a maximum. In this letter, using the first data from 



The discovery that approximately half the energy ever radiated the Herschel Multi-tiered Extragalactic Survey (HerMES; Oliver 

J-j ; by galaxies is received on Earth in the far-infrared waveband et ^ in }> we inves tigate the evolution implied by the ex- 

^N , (Puget et al. 1996; Fixsen et al. 1998) implies that galaxies must istence of the FIRB by measuring the evolution of the galaxy 

IQ show strong evolution that is hidden from optical telescopes luminosity function at 2 50 pm. We everywhere assume a stan- 

O (Gispert, Lagache & Puget 2000). A decade ago, deep surveys dard concordance cosmology: Q M = 0.28, Q A = 0.72, H„ = 

^^ with the SCUBA camera on the James Clerk Maxwell Telescope 7? km s ' Mdc ' 
i~ I ' resolved much of the far-infrared background (FIRB) at 850 pm 
£>. | into individual sources (Barger et al. 1998; Hughes et al. 1998). 

These sources are mostly extremely luminous dust-enshrouded 2 The data 
galaxies at z > 2 (Chapman et al. 2005) with an average im- 
plied star-formation rate (if the ultimate source of the energy is The images that we analyse in this letter were taken at 250 pm 
star formation) of ^400 M year" 1 (Coppin et al. 2006), much with the SPIRE instrument on Herschel, whose in-orbit perfor- 
greater than the star-formation rates in galaxies like our own. mance and scientific capabilities are described in Griffin et al. 
However, the energy density in the FIRB at 850 pm is ^30 (2010). The calibration methods and accuracy of SPIRE are de- 
times less than at 200 pm where the FIRB is at a maximum, scribed by Swinyward et al. (2010). The three images consist of 
and both the spectral shape of the FIRB and statistical 'stack- a shallow image of the Lockman Hole (LH) and deep images 
ing' analyses (Dole et al. 2006; Pascale et al. 2009) imply that of the northern field of the Great Observatories Origins Deep 
much of the FIRB is actually produced by sources at lower red- Survey (GOODS-North) and of a field within the Lockman Hole 
shift (Gispert et al. 2000; Dole et al. 2006; Pascale et al. 2009). (LH-North). For the latter two images the dominant source of 
The launch of the Herschel Space Observatory (Pilbratt et al. noise is confusion due to numerous faint sources, which is 5.8 
2010) in May 2009 has given us the opportunity to resolve a mJy beairT 1 at 250 yum (Nguyen et al. 2010). 

To measure robust fluxes for sources close to the confusion 
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* Herschel is an ESA space observatory with science instruments level, we have developed a source-extraction technique that is 
provided by European-led Principal Investigator consortia and with im- based on the assumption that sources detected at 250 pm with 
portant participation from NASA. SPIRE will also be detected in deep observations with Spitzer 
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Fig. 1. Fraction of sources recovered by the source-extraction 
method as a function of flux density for GOODS-North (dashed), 
LH-North (dot-dashed) and the LH (dotted). The vertical lines 
show the flux limits we have used for the samples from the three 
fields (Table 1). 

Table 1. Field Statistics 



Field 


GOODS-North 


LH-North 


LH 


Area (deg-) 


0.031 


0.371 


5.17 


Flux limit (mjy) 


10 


20 


48 


No. sources 


83 


276 


551 


Spectroscopic z 


65 


75 


140 


Photometric z 


18 


161 


322 


No z (%) 





14 


16 



at 24 pm. This assumption was suggested by the recent studies 
that concluded that galaxies detected in deep 24-jum surveys with 
Spitzer produce most of the FIRB at 160 pm (Dole et al. 2006) 
and at 250, 350 and 500 pm (Marsden et al. 2009). This approach 
reduces the effective confusion noise by resolving some of the 
confusing background into the sources detected at 24 pm. Our 
'cross-ID method' starts from a list of the 24-yum sources found 
in the field covered by the 250-jt/m image. By using a matrix 
inversion technique, we then find the 250-pm fluxes at these po- 
sitions that provide the best fit to the 250-/vm image (Roseboom, 
in prep). A problem with methods like this is that there can be 
a large number of degenerate solutions when the surface-density 
of sources in the input catalogue is large, as it is with the deep 
Spitzer 24-yum catalogues. We have addressed this problem by 
iteratively reducing the number of sources in the input catalogue 
in order to find the list of input sources that provides the best 
fit to the 250-/vm image. Full details of the method, its valida- 
tion with simulations, and a comparison of the results with other 
methods of source extraction are given in Roseboom et al. (in 
prep). 

Roseboom et al. (in prep) have assessed the completeness of 
the cross-ID method by inserting artificial sources onto the real 
images and then determining the fraction of these sources that 
are detected by the source-extraction method. Fig. 1 shows the 
results for the three fields as a function of flux density. Using 
these curves to choose the flux limits, we have used the cross-ID 
method to extract samples of sources from regions within each 
field for which the optical/IR data is particularly good (Table 
1). We accepted a lower completeness level at the flux limit 
for GOODS-North (^50%) in order to get the maximum range 
of luminosity at each redshift. When calculating the luminosity 
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Fig. 2. The histogram shows the distribution of 24-fim flux den- 
sity for the 5015 24-/im sources in the Spitzer catalogue used in 
the source extraction for LH-North. The red line shows the dis- 
tribution of 24-pm flux density for the sample of 276 250-yi/m 
sources found using our source-extraction method. 



function (§3), we used the curves in Fig. 1 to correct for incom- 
pleteness. All the sources above these flux limits are detected at 

> 5cr. 

Fig. 1 does not account for any objects that are missing be- 
cause their 24-pm flux density falls below the limit of the Spitzer 
catalogue used in the cross-ID method. There are two arguments 
that suggest this is not a serious problem. First, Roseboom et al. 
show that the number-density of 250-yum sources in the cross- 
ID catalogues agree well with the source counts determined by 
Oliver et al. (2010) down to ^25 mjy in LH-North and ^40 mjy 
in LH. Second, they compare the results of using a shallow and 
a deep Spitzer catalogue in several fields and conclude that us- 
ing the deeper catalogue does not produce a large increase in the 
number of 250-fim sources found by the cross-ID method. Fig. 2 
illustrates this clearly for the LH-North field. The figure shows 
a histogram of the 24-pm flux densities of the sources in the in- 
put Spitzer catalogue, whereas the red line shows the 24-pm flux 
densities of the objects found by the cross-ID method. It is im- 
portant to note that the 24-pm flux densities of the sources in 
the input catalogue are not used as information in the cross-ID 
method, so the very different distributions for the input catalogue 
and for the objects found by the cross-ID method is evidence that 
the method is working well. The small number of objects with 
faint 24-pm flux densities found by the cross-ID method, de- 
spite the very large number of faint 24-pm sources in the input 
catalogue, strongly suggests that we are not missing objects that 
have too low 24-pm flux densities to be included in the Spitzer 
catalogues. 

We have used the optical and near-infrared images for these 
fields to find the counterparts to the 24-pm sources. Table 1 
shows we have redshifts for all of the GOODS-North sources, 
either a spectroscopic redshift from the catalogue of Barger et al. 
(2008) or, for 22% of the sources, a photometric redshift, which 
we have estimated from the available multi-wavelength images 
(in typically nine optical and near-infrared bands - Raymond 
et al., in prep). The situation for our other deep field is a little 
less satisfactory, although we still have redshifts for 86% of the 
sources, a mixture of spectroscopic redshifts and photometric 
redshifts, mostly taken from an unpublished catalogue of Owen 
and collaborators, which was produced from images in eight op- 
tical and near-infrared bands. The data for LH is described in 
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Fig. 3. Redshift histograms for GOODS-North (continuous) and 
LH-North (dashed). The values in the histogram for LH-North 
have been scaled to the number of sources found in GOODS- 
North. 
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Fig. 4. Luminosity function at 250 //m in six redshift slices: 
0.0 < z < 0.2— red; 0.2 < z < 0.4— green; 0.4 < z < 0.8— 
dark blue; 0.8 < z < 1.2 — light blue; 1.2 < z < 1.6 — purple; 
1.6 < z < 2.0 — yellow. The black diamonds show the esti- 
mate of the local 250-/itn luminosity function from Vaccari et 
al. (2010). 



Rowan-Robinson et al. (2008) and Vaccari et al. (in prep). Fig. 3 
shows the redshift distributions for our two deep samples. The 
differences between the two can probably be explained by cos- 
mic variance as the result of the small area of the GOODS-North 
sample. The similarity of the distributions at high redshift sug- 
gests that the sources without redshifts in LH-North are not pref- 
erentially at high redshift. 



3. Estimates of the luminosity function 

We have estimated the rest-frame 250-/im luminosity function 
in six redshift intervals: 0.0 < z < 0.2, 0.2 < z < 0.4, 0.4 < 
z < 0.8, 0.8 < z < 1.2, 1.2 < z < 1.6, 1.6 < z < 2.0. Because 
the optical/IR data is less sensitive for this field than the other 



two, we have only used the LH sample for the two low-redshift 
intervals. We have not estimated the luminosity function at z > 2 
because of the low number of sources above this redshift. 

At this early stage of the HerMES project, we do not have 
large numbers of measurements of the spectral energy distribu- 
tions (SED) of individual galaxies. We have therefore made the 
assumption that all galaxies have the same SED, which we have 
taken to be a grey body with a dust-emissivity index of 1 .5 and 
a temperature of 26 K, the average of the SEDs found for the 
galaxies detected at the same wavelength by the Balloon-borne 
Large-Aperture Sub-millimeter Telescope (BLAST; Dye et al. 
2009). If the typical temperature of the dust in a SPIRE galaxy 
does change systematically with redshift, this will have little ef- 
fect on the shape of a galaxy SED at A > 100 /mi and thus on 
the luminosity functions at z < 1 but it might have a significant 
effect on the luminosity functions at higher redshift. 

We have estimated the luminosity function using the esti- 
mator of Page and Carrera (2000), which is not only the ideal 
estimator in general but is particularly well-suited to situations 
where source confusion is an issue (Eales et al. 2009). For each 
source, we used the curves in Fig. 1 to estimate a correction fac- 
tor, Ci, which is the reciprocal of the completeness at the flux 
density of the source. Our estimate of the value of the lumi- 
nosity function in a bin of the luminosity-redshift plane is then 
2 Cj/V, where the sum is over the galaxies in this bin and V is 
the accessible volume in this bin averaged over the width of the 
bin in luminosity. We have combined the data for the different 
HerMES fields by adopting the 'coherent volume' approach of 
Avni and Bahcall (1980). Our estimate of the luminosity func- 
tion in a bin is then 2 Q/ Z V, where the sum in the numerator 
is now over all galaxies in the three fields in that bin and the 
sum in the denominator is the sum of the accessible volumes for 
the three fields. The fractional error on the luminosity function is 
then yN/N, in which N is the uncorrected number of galaxies in 
that bin in the three fields. This error does not include the effect 
of cosmic variance, which we will consider in a later paper. 

Fig. 4 shows our estimates of the luminosity function in the 
six redshift intervals and the HerMES estimate of the local lumi- 
nosity function (Vaccari et al. 2010), which was estimated over 
two fields covering 14.7 deg 2 (including our LH area) and agrees 
reassuringly with our estimate in the lowest redshift slice. The 
luminosity function shows steady evolution out to z = 1, in the 
sense that the space-density of the most luminous sources grad- 
ually increases with redshift. This agrees well with the evolution 
in the luminosity function determined from the BLAST results 
(Eales et al. 2009). There is evidence in Fig. 4 for at most weak 
evolution at z > 1 . 

There are a number of possible systematic errors that could 
affect our results. One of the nice properties of the cross-ID 
method is that simulations suggest it is relatively immune to the 
effect of flux boosting due to source confusion and Eddington 
bias (Roseboom et al., in prep). A second possible problem is 
that we do not have redshifts for all the sources in LH and LH- 
North. If the incomplete redshift information affects each red- 
shift bin equally, the luminosity functions will all be slightly 
too low. However, if the incompleteness preferentially affects the 
high-redshift bins, we will have under-estimated the size of the 
evolution. Another potential problem is the large number of pho- 
tometric redshifts in our analysis. We have used a Monte-Carlo 
simulation to generate artificial samples based on a no-evolution 
assumption, giving each source a photometric redshift with an 
accuracy equal to our least accurate redshifts, and then used the 
artificial samples to estimate the luminosity function in the six 
redshift slices. This analysis confirms that the large number of 
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photometric redshifts does not generate spurious strong evolu- 
tion. Finally, even if the dust in galaxies does not get system- 
atically hotter with redshift, there will be a selection effect in 
our results, since at high redshift we are sampling a lower rest- 
frame wavelength, and the galaxies found in samples selected 
at shorter wavelengths tend to contain hotter dust. The effect of 
this selection effect is that we may have overestimated slightly 
the luminosity functions at z > 1 . 

4. Discussion 

On the assumption that the star-formation rate is proportional 
to the rest-frame 250-yum luminosity, the lack of evidence for 
any strong evolution at z > 1 is consistent with investigations 
that have concluded that the overall star-formation rate per unit 
comoving volume was approximately constant at z > 1 (Steidel 
et al. 1999; Gispert et al. 2000; Hopkins 2004). We estimate from 
the GOODS-North catalogue that we have resolved ^20% of the 
FIRB at 250 /im. Fig. 3 shows that the sources making up this top 
20% of the FIRB are at moderate redshift (z ~ 1). Although we 
cannot say anything directly about the sources responsible for 
the missing 80% of the FIRB, the stacking analysis of Pascale 
et al. (2009) from BLAST data suggests that the remainder of 
the FIRB at this wavelength is also from sources at moderate 
redshift. 

A revealing way to look at these early HerMES results is 
to use the relationship between the production of metals and 
the background radiation associated with this metal produc- 
tion (Peacock 1993). Using a value for the integrated FIRB of 
14 nW sr~' (Fixsen et al. 1998), we obtain a relationship between 
the mass of metals and the fraction, e, of the FIRB produced in 
a particular redshift interval: 

M melals = 9 x 10 6 e(l + z) M Q Mpc-\ 

Fig. 5 shows the metals produced as a function of redshift on 

the assumption that the redshift distributions of GOODS-North 
and LH-North are representative of the FIRB as a whole. The 
figure suggests that most of the metals (and therefore most of 
the massive stars) formed at a moderate redshift (z < 1.4). This 
conclusion is consistent with the overall star-formation rate in 
the universe being constant at z > 1 , and the apparent decline of 
metal production at high redshift is simply the consequence of 
the relationship between cosmic time and redshift. However, the 
decline in the figure at high redshifts is relatively small, and since 
we are making the large assumptions that the redshift distribu- 
tions of our samples are characteristic of the FIRB as a whole, it 
is possible that observations that resolve more of the FIRB and at 
several wavelengths rather than a single wavelength will modify 
this conclusion. 

The excellent images that exist from the Hubble Space 
Telescope for the GOODS-North Herschel sources mean that we 
can examine the nature of the galaxies responsible for this metal 
production. We have used the z-band images, since at z =* 1 
this band corresponds in the rest-frame to the B-band. Of the 
26 galaxies in the redshift interval 0.8 < z < 1.2, 12 show 
clear signs of spiral structure, and there may well be more if 
some of the galaxies have spiral arms that are below the surface- 
brightness detection threshold of the images. Thus the galaxies 
found in the first HerMES images appear to be quite different 
from the major mergers found in the SCUBA surveys (Ivison et 
al. 2000; Tacconi et al. 2008). Studies of galaxy evolution based 
on Spitzer 24-yt/m surveys have also found that the Spitzer galax- 
ies at z =* 1 are often spirals rather than major mergers (Elbaz et 




Redshift 

Fig. 5. An estimate of the total metals formed per cubic Mpc 
in different redshift intervals using the method described in the 
text. We have made the estimates from the redshift distributions 
of GOODS-North (red) and LH-North (green) on the assumption 
that these redshift distributions represent the FIRB as a whole. 



al. 2007) but the crucial advance that Herschel has made possible 
is that we can determine the cosmological importance of these 
objects, in terms of the total masses of stars and metals formed 
in them. 
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